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ABSTRACT 

We have developed a high resolution combined physical and chemical model of a protoplanetary 
disk surrounding a typical T Tauri star. Our aims were to use our model to calculate the chemical 
structure of disks on small scales (sub-milli-arcsecond in the inner disk for objects at the distance of 
Taurus, ^ 140 pc) to investigate the various chemical processes thought to be important in disks and 
to determine potential molecular tracers of each process. Our gas-phase network was extracted from 
the UMIST Database for Astrochemistry to which we added gas-grain interactions including freeze 
out and thermal and non-thermal desorption (cosmic-ray induced desorption, photodesorption and 
X-ray desorption) and a grain-surface network. We find that cosmic-ray induced desorption has the 
least effect on our disk chemical structure while photodesorption has a significant effect, enhancing 
the abundances of most gas-phase molecules throughout the disk and affecting the abundances and 
distribution of HCN, CN and CS, in particular. In the outer disk, we also see enhancements in the 
abundances of II2O and CO2. X-ray desorption is a potentially powerful mechanism in disks, acting to 
homogenise the fractional abundances of gas-phase species across the depth and increasing the column 
densities of most molecules although there remain significant uncertainties in the rates adopted for this 
process. The addition of grain-surface chemistry enhances the fractional abundances of several small 
complex organic molecules including CH3OH, HCOOCH3 and CH3OCH3 to potentially observable 
values (i.e. a fractional abundance of > 10^^^). 

Subject headings: astrochemistry — planetary systems: protoplanetary disks — stars: formation — 
solar system: formation — ISM: molecules 



1. INTRODUCTION 

Protoplanetary disks are crucial objects in low-mass 
star formation, possessing three vital functions: they (i) 
aid the dissipation of angular momentum away from the 
young stellar system, (ii) allow the efficient accretion of 
matter onto the young star and (iii) contain all material, 
dust and gas, which may end up in a planetary system 
orbiting the main-sequence star. 

In this work, we investigate the chemistry and molec- 
ular composition of a protoplanetary disk surrounding 
a young star which will evolve into a main-sequence 
star resembling our Sun. At the low temperatures 
encountered in many astrophysical regions (~ 10 K 
to ~ 100 K), molecules are readily excited into higher 
rotational energy states and subsequently emit radiation 
at (sub)millimeter wavelengths. Early observations 
of T Tauri stars at these wavelengths revealed the 
presence of m olecular mater i al in a flattened disk-like 
structure (e.g. iDutrev et al.l ()1994D) and in Keplerian 
rotation about the parent star (e.g.' Guilloteau &: DutrevI 
Since then, molecular rotational line emission 
originating from a disk has been obser yed in several 
T Ta u ri sy stems inclu ding TW H y drae (^Kastne r et"aL 
19971: Jwi Zadelhoff e t"al1 120011: Ivan Dishoeck et al. 
2003; 'Ceccarelh et alj 120041: iThi et al.' 2004: Oi et al. 
2004 , 2006, 20081 D M Ta u ri (iDutrev et al 19.s5 
Ceccarelli et al.l l200l 120051 : IGuiUoteau et all 12006 



IDutrev et al.l [20071: iPietu et al.l [20071) and Lk Ca 15 



van Zad elhoff et al.l ffoOlt lAikawa et al.l [2 003: Oi et al 



iThi et all 120041: IDutrev et al.l 12007c iPietu et al 
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12007!). Most species detected are small simple molecules, 
molecular ions and radicals such as CO, HCO"*", CN, 
HCN, CS, C2H and N2H+, along with several associated 
isotopologues (e.g. ^^CO, C^^O, DCO+, H"CN, H2D+ 
and C'^^S). The most complex species observed to date 
is the small organ i c molecule, formald e hyde, H2CO 
(IDutrev et al.l 119971: lAikawa et all 120031: IDutrev et akl 
12007!) with methanol C H3OH, thus far e luding detection 
(e.g. lThi et all (|200l ). iCeccarelh et"all ([2005) report a 
detection of deuterated water, HDO, in the disk of DM 
Tau, although this re sult has since been disputed by 
'Guillo teau et al.l (|2006[ ). 

Infra-red emission has also been observed originating 
from disks embedded in young stellar objects and aris- 
ing from vibrational transitions in gas-phase molecules 
capable of survival in the warmest regions (> 350 K). 
Thus, infra-red emission probes not only a different phys- 
ical region of the disk to that probed by (sub)mm emis- 
sion, but also uses different molecules as tracers, hence 
providing complimentary chemical information. The 
molecules detected thus far at infra-re d wavelengths are 



CO, HCN, OH, H^O , CO^ 
Lahuis et al. 2006j_ 
2008t iPascucci et al 



: and C^H9 (ICarr et al.l[200l 
ICarr fc Naiit^ 120081 : ISalvk et al.l 
2009D with an upper limit deter- 
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Observations of molecular line emission from disks, to 
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date, have been hampered by the small angular size of 
these objects on the sky and the limitations of exist- 
ing facilities, explaining why the species detected are 
those which are abundant and possess relatively simple 
rotational energy spectra (e.g. CO). Single-dish facilities 
which operate at (sub)mm wavelengths such as the 15 m 
James Clerk Maxwell Telescope (JCMT) in Hawaii and 
the IRAM 30 m telescope in Spain, have been predomi- 
nantly used in the detections of the molecular species in 
the T Tauri systems listed. With beam-sizes much larger 
than the typical source size, usually a single molecular 
line profile is generated characterising emission from the 
entire disk. In order to spatially resolve the emission and 
hence trace the radial and vertical physical and chemical 
structure , interferometry must be employed and indeed, 
iQi et al.l J2008) report spatially resolved emission arising 
from molecular rotational transitions in the disk of TW 
Hya using the Sub-Millimeter Array (SMA). 

The discipline of (sub)mm astronomy is sched- 
uled for a revolutionary transformation with the 
first light of the Atacama Large Millimeter Ar- 
ray (ALMA) in Chile expected in 2012 (see 
http://www.almaobservatory.org). ALMA, with 
its 50 12 m telescopes and fully variable configuration, 
will have the spatial resolution necessary to observe 
molecular line emission from protoplanetary disks on 
sub-milli-arcsecond scales and enable the tracing of the 
molecular content of disks to within « 0.1 AU of the 
parent star at its highest operational frequencies. It 
is anticipated that the sensitivity and high spectral 
resolution of ALMA will lead to the potentially over- 
whelming detection of many further molecular species, 
including complex organic molecules considered the 
building blocks of life, in many astrophysical sources 
including protoplanetary disks. 

Motivated by the impending completion of ALMA, we 
have constructed a high resolution combined chemical 
and physical model of a protoplanetary disk surrounding 
a typical T Tauri star using as comprehensive a chem- 
ical network as computationally possible. In the work 
presented here, our objectives were (i) to calculate the 
chemical structure of protoplanetary disks on small (sub- 
milli-arcsecond in the inner disk) scales, (ii) to investigate 
the infiuence of various chemical processes, such as non- 
thermal desorption and grain-surface chemistry, thought 
to be important in disks, and (iii) to subsequently de- 
termine potential molecular tracers of each process. We 
also used our model to (i) compute molecular line emis- 
sion profiles for rotational transitions which have been 
observed in disks using existing facilities, (ii) compare 
our modelled line profiles and intensities with existing 
observations and (iii) produce molecular line emission 
maps at the expected spatial resolution of ALMA for 
disks at various distances and inclinations. This second 
study and corresponding set of results will be covered 
in a subsequent paper (Walsh et al. in preparation). 
Our study also aims to help answer some fundamental 
questions concerning the evolution of stars, planets and 
ultimately, life. Is it possible for primordial (possibly or- 
ganic) material created in a young star's protoplanetary 
disk to survive the assimilation into planets and other 
planetary system objects? Is our solar system's chemical 
and thus, planetary composition unique? How intrinsi- 
cally linked are star formation, planet formation and the 



origin of life in the universe? These questions are ev- 
ermore important as we move into the era of exoplanet 
research and the hunt for planets and the signatures of 
life in external stellar systems. 

In Section[2]we describe the theoretical foundation and 
generation of the physical model used to characterise 
our protoplanetary disk (Section 12. ip and the chemi- 
cal network we have collated and used in our calcula- 
tion of the disk chemical evolution (Section 12. 2p includ- 
ing gas-phase chemistry (Section l2.2.ip . photochemistry 
(Section l2.2.2[) . gas-grain interactions (Section l2.2.3p and 
grain-surface chemistry (Section I2.2.4p . The results of 
our chemical evolution calculations are covered in Sec- 
tion [3] where we discuss the chemical structure and strat- 
ification in the disk (Section 13. ip . the effects of our in- 
cluded chemical processes (Sections 13. 2l and l3.3p . the disk 
ionisation fraction (Section 13.41) and the radial molecu- 
lar column densities (Section 13. 5p . We briefly discuss 
our work in relation to similar projects by other research 
groups in Section 13.61 and finally, in Section |4l we sum- 
marise our work and outline our main conclusions and 
further work we intend to undertake. 

2. PROTOPLANETARY DISK MODEL 
2.1. Physical Model 

The physical mod el of a protoplanetary di sk we use 
in this work is from INomura fc Millail (I2005D with the 
addition of X-ray heating as described in INomura et al.l 
([2007). They self-consistently modelled the density and 
temperature profiles of gas and dust in a protoplane- 
tary disk accounting for UV and X-ray irradiation by 
the central star and subsequently computed molecular 
hydrogen line emission at ultraviolet and infrared wave- 
lengths. Here, we have used this model to compute the 
chemical structure of a protoplanetary disk with the ulti- 
mate aim to expand on their work by calculating molecu- 
lar line emission from disks at (sub)mm wavelengths. In 
the remainder of this section, we give a brief overview of 
our physical model and we refer readers to the original 
papers for the mathematical and computational details. 

We consider an axisymmetric disk surrounding a typ- 
ical T Tauri star with mass, Af* = 0.5 Mq, ra- 
dius, = 2 Rp; and te mperature, = 4000 K 
(jKenvon fc HartmannI 119 95^. The density and temper- 
ature distributions are determined through iteratively 
solving the equations for hydrostatic equilibrium in the 
vertical direction and the local thermal balance between 
the heating and cooling of the gas. The theoretical foun- 
dation of this model comes from the stand ard a ccretion 
disk m odel of iLvnden-Bell fc Pringld (|1974l ) and iPringli 
(|1981| ) which defines a surface density distribution for 
the disk given the parent star's mass and radius and a 
disk mass accretion rate, M. The kinematic viscosity 
in the disk is parameteris ed according to the work of 
iShakura Sz SunvaevI (|1973f ). the so-called a-prescription. 
We adopt a viscous parameter, a — 0.01 and a mass 
accretion rate, M — 10~* Mq yr~^. 

The heating mechanisms included are grain photo- 
electric heating by far-ultraviolet photons and X-ray 
heating due to hydrogen ionisation by X-ray photons 
with cooling via gas-grain collisions and line transitions. 
We use a model spectrum created by fitting the observed 
XMM-Newton X-ray spectrum of the classical T Tauri 
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star, TW Hya (e.g. iKastner et all (|2002D ) with a two- 
temper ature thin thermal pla sma model (MEKAL model 
see e.g. iLiedahl et al.l (|1995[ )). The X-ray luminosity is 
Lx ~ 10'^° erg s~^ and the resulting X-r ay spectrum is 
given in Figure 1 of iNomura et al.l (|2007D . 

The UV radiation field in disks has two sources, the 
star and the interstellar medium. In this disk model, 
the radiation field due to the T Tauri star has three 
components: black-body emission at the star's efi'ective 
temperature, optically thin hydrogenic bremsstrahlung 
emission and strong Lyman-a line emission. All compo- 
nents are necessary to accurately model the excess UV 
emission observed towards classical T Tauri stars thought 
to arise from an accretion shock as disk material im - 
pinges upon the stellar surface (e.g. lValenti et al.l (|2000l ): 
[johns-K ruU et al.l (2000)). The total FUV luminosity in 
our model is Ljjv ~ 10'^^ erg with the calculation of 
the radiatioi i field in the disk describ ed in detail in Ap- 
pendix C of INomura fc Millail ()2005D and the resulting 
spectrum shown in Figure C.l in that paper. 

We assume the dust and gas in the disk is well-mixed 
and adopt the dust-size distribution model which repro- 
duces the observational extinction curve of dense clouds 
(jWeingartner fc Draind l2001| ) . The calculation of the 
dust opacity in th e disk is as described in Appendix D of 
INomura fc Millail (2005) with the resulting monochro- 
matic absorption coefficient shown in Figure D.l. We 
note here that this is an over-simplification of the treat- 
ment of the dust-size distribution in protoplanetary disks 
and we are currently working on improving our model by 
adding in the effects of dust-grain settling and coagula- 
tion. 

In Figure [T] we display the resulting number density 
(cm~^), gas temperature (K) and dust temperature (K) 
as a function of disk radius and height (top, middle 
and bottom rows, respectively). To illustrate the ex- 
treme vertical gradients in the physical conditions at 
small radii, we display the density and temperature both 
within 10 AU (left panels) and 305 AU (right panels). 
We describe the physical structure of our disk in the Ap- 
pendix. 

2.2. Chemical Model 

The structure of the disk described in the preceding 
section leads to a multitude of different physical regimes 
and as such, we need to account for every chemical pro- 
cess which may occur. The axisymmetric structure re- 
sults in a cold, dense midplane where even the most 
volatile molecules are expected to freeze out onto dust 
grains creating an icy mantle and depleting the gas of 
molecules. Moving in the vertical direction, the den- 
sity decreases and the temperature increases driving the 
evaporation of molecules from grain surfaces and stim- 
ulating a rich gas-phase chemistry resulting in further 
molecular synthesis. Further towards the surface, the 
radiation fields increase in strength dissociating and ion- 
ising molecules into constituent radicals, atoms and ions. 
A similar stratification is expected in the radial direction 
as the temperature and density in the disk midplane both 
increase with decreasing distance from the star. When 
the midplane dust temperature reaches a value higher 
than the desorption temperature of a particular molecule, 
it is returned to the gas phase. This point is known as 



the snow line and can occur at a unique radius for each 
molecule. At small radii, due to the high densities found 
in the midplane, there is a significant column density of 
material shielding this region from the intense UV and 
X-ray fields of the star such that molecules are expected 
to survive in the midplane at radii within ~ 0.1 AU. 

In order to investigate the chemical structure thor- 
oughly, we used a large gas-phase network supplemented 
with gas-grain interactions, including freeze out and ther- 
mal desorption. We considered various non-thermal des- 
orption mechanisms, namely, cosmic-ray induced desorp- 
tion, photodesorption and X-ray desorption. To probe 
the efficacy of molecular synthesis on grain-surfaces we 
also added a large grain-surface reaction network. 

2.2.1. Gas-Phase Chemistry 

Our gas-phase chemistry is extracted from the latest 
release of the 'dipole-enhanced' version of the UMIST 
Database for Astrochemistry ( http : // www . udf a . net ) , 
henceforth referred to as 'Rate06' (Woo dall et al.l [2007). 
We include almost the entire Rate06 gas-phase network 
removing only those species (and thus reactions) which 
contain fluorine, F, and phosphorus, P, in order to re- 
duce computation time. We deemed the loss of F- and 
P-containing species to have a minimal impact on the re- 
maining chemistry. Our gas-phase network thus consists 
of 4336 reactions involving 378 species composed of the 
elements H, He, C, N, O, Na, Mg, Si, S, CI and Fe. The 
initial elemental fractional abundances (relative to total 
H nuclei density) we use ar e the set of oxygen-ri ch low- 
metallicity a bundances from[Graedel et al.l (|1982l ). listed 
in Table 8 of lWoodall et~all ^Wf). We flnd by 10^ years, 
the typical age of protoplanetary disks, the chemistry has 
forgotten its origins, justifying our use of initial elemen- 
tal abundances. We intend in future models to calculate 
the chemical evolution of a parcel of gas as it follows 
a streamline in the accretion flow in which case the in- 
put abundances should reflect the molecular make-up of 
the ambient cloud material. Our model grid has over 
12,000 grid points in 129 logarithmically spaced radial 
steps from 0.04 AU to 305 AU. 

2.2.2. Photochemistry 

In the models presented here, we have approximated 
our photoreaction rates at each point in the disk, 
fcP'*(r, z), by scaling the rates from Rate06 (which as- 
sume the interstellar UV field) using the wavelength in- 
tegrated UV flux calculated at each point, Gpuvif, z) — 

-CiTl"^ Gi?t/y (A, r, z) dA. Hence, the rate for a particular 
photoreaction at each (r, z) is given by 

F'^^^fco B-\ (1) 

where Go is the interstellar UV flux and fco is the rate 
expected in the interstellar medium. 

2.2.3. Gas-Gram Interactions 

Gas-grain interactions are important in large areas 
of protoplanetary disks as the dust temperature can 
reach values lower than the freeze-out temperatures of 
molecules. If the freeze out of gas-phase species is al- 
lowed, then the evaporation of molecules from dust grains 
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Figure 1. Number density (top), gas temperature (middle) and dust temperature (bottom) as a function of disk radius and height up to 
maximum radii of r = 10 AU (left) and 305 AU (right). 



must also be included. In this work, we consider both the 
thermal and non-thermal desorption of molecules from 
dust grains. For the thermal desorption of a particular 
molecule to occur, the dust-grain temperature must ex- 
ceed the freeze-out temperature of that molecule. Non- 
thermal desorption requires an input of energy from an 
external source and is thus independent of dust-grain 
temperature. As protoplanetary disks are irradiated by 
UV and X-ray photons from the central star as well as 
UV photons and cosmic-rays originating from the in- 
terstellar medium, the non-thermal desorption mecha- 
nisms we investigate are cosmic-ray induced desorption, 
photodesorption, and X-ray desorption. Our gas-phase 
chemical network has thus been supplemented with an 
additional 1154 gas-grain interactions involving 149 sur- 
face species. 



kf, of species 



The accretion rate (or freeze-out rate) 
i onto dust-g rain surfaces is treate d using the standard 
prescription (|Hasegawa et al.lll992f ). 



(2) 



where Si is the sticking coefficient, here assumed to 
equal unity for all species, ad = tto^ is the geometri- 
cal cross section of a dust grain with radius, a, < Vi >= 
(kBT/rrii)^^'^ is the thermal velocity of species j at a tem- 
perature, T and with mass, rrii, ks is Boltzmann's con- 
stant, and Ud is the number density of dust grains. 

The thermal desorption rate, A;f, of species i is de- 
pendent on dust-grain temperature, Td (|Hasegawa et al.l 



Chemical Processes in Protoplanetary Disks 



5 



[T99l) . and is given by 

kf = voii) exp 



-Edit) 
Td 



(3) 



where Ed[i) is the binding energy of species i to the dust 
grain in units of K. The characteristic vibrational fre- 
quency of each adsorbed species in its potential well, 
vi^iji), is represented by a harmonic oscillator relation 
(jHasegawa et al.lll992[ ). 



'2nsEd{i) 



(4) 



where, here, Ed{i) is in units of erg and Us ~ 1.5 x 
10^^ cm^^ is the number density of surface sites on 
each dust grain. The binding energies, Ed-, for sev- 
er al important molecules (mai nly following those co llated 
bv IHasegawa et all (|1992[ ) and lWillacv et all (|1998 )) arc 
listed in Table [TJ We intend to conduct a review of our 
set of desorption energies in light of more recent experi- 
mental results e.g. a recent investigation in to the desorp- 
tion of methanol bv lBrown fc Bolinal (|2007f) determined a 
binding energy of w 5000 K, as opposed to the theoretical 
value of 2140 K used here (see Table [T]). This binding en- 
ergy was determined for pure methanol ice as opposed to 
methanol adsorbed onto, or mixed with, water ice. We 
find throughout our model that the ratio of methanol 
to water ice is less than 1%. Similar experiments for 
methanol adsorbed onto water ice (Brown, private com- 
munication) suggest that the binding energy of methanol 
in this complex is comparable with that determined for 
pure methanol but due to overlapping desorption fea- 
tures the results are d ifficult to analyse. Recent work by 
iBottinelli etldl (|2010l ) comparing laboratory data with 
observations of methanol in young stellar objects (YSOs) 
suggests that methanol ice in these environments likely 
exists as pure ice or mixed with CO and/or CO2 ice which 
is consistent with its formation via hydrogenation of CO 
on dust grains. Considering the latter molecules are non- 
polar, it is possible that the binding energy of methanol 
in astrophysical ices is lower than that determined in the 
laboratory experiments. Increasing the binding energy of 
methanol to a value of ~ 5000 K will increase the desorp- 
tion temperature from f« 30 K to 40 K to ~ 100 K. This 
will push the 'snow line' for methanol closer to the star 
but should have little effect on the outer disk methanol 
abundances where the dust temperature is < 30 K. We 
expect non-thermal desorption to dominate over thermal 
desorption in the upper layers of the disk. 

To calculate the cosmic-ray induced des orption rate 
for ea ch spe cies, fcf'^, we use the method of ILeger et al.l 
(fl985l ) and IHasegawa fc HerbstI (fl993l) . They assume 
that dust grains with a radius of 0.1 /zm are impulsively 
heated by the impact of relativistic Fe nuclei with ener- 
gies of 20 to 70 MeV nucleon"^ which deposit, on aver- 
age, an energy of 0.4 MeV into each dust grain. Assum- 
ing that the majority of molecules desorb around 70 K, 
the cosmic-ray induced desorption rate can be approxi- 
mated by 



/(70 K)fcf (70 K) s" 



(5) 



where kf (70 K) is the the thermal desorption energy of 
species i at a temperature of 70 K, calculated using Equa- 



Table 1 

Molecular Bindine 



Energies 



Molecule Binding Energy (K) Reference 



CO 

N2 

HCN 

CO2 

H2O 

NH3 

CH4 

C2H2 

H2CO 

CH3OH 



960 
710 
4170 
2690 
4820 
3080 
1080 
2400 
1760 
2140 
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tion The parameter, /(70 K), is the fraction of time 
spent by grains in the vicinity of 70 K and can loosely 
defined as the ratio of the desorption cooling time (sa 
10"^ s~^) to the time interval between successive heat- 
ings to 70 K (3.16 X 10 " s) so that f (70 K) » 3.16x lO'^^ 
(for further details see IHasegawa fc Herbst] (|1993D ). 

Note that the method of calculating the cosmic-ray in- 
duced desorption rates is species dependent and a func- 
tion of surface binding energy. In contrast, the photodes- 
orption rates are indiscriminate, based on the experi- 
mental results of iWestlev et al.l (jl995') and lOberg et al.l 
(2007). Their results suggest each photon absorbed 
by the grain mantle returns a particular number of 
molecules, independent of binding energy, to the gas 
phase so that the desorption rate of each species varies 
according to its fractional abundance on dust-grain sur- 
faces. The overall phot odesorption rate is calcu lated, 
similar t o the work of iWillacv fc Langeil ()2000[ ) and 
iWillacvl (l2007l) . using 

kP'' ^ FuvYuv^dXd s-\ (6) 

where, Fijv is the UV radiative flux in units of pho- 
tons cm^^ s^"'^, Yuv is the experimentally determined 
photodesorption yield in units of molecules photon"^, 
(Jd is the geometrical dust-grain cross section in cm^ and 
Xd is the fractional abundance of dust grains. Note that 
the attenuation of UV radiation is accounted for in our 
calculation of Fjjv- Hence, the photodesorption rate for 
a specific species, kf^, is calculated using 



(7) 



where, feP"* is given by Equation ([6]), nf is the number 
density of species i frozen out onto grain surfaces and 
nf^j is the total number density of grain-surface species. 

More recent experiments bv lOberg et al.l ([200930) sug- 
gest that photodesorption rates are also dependent on the 
depth of the ice layer on grain surfaces with the molec- 
ular yield also dependent on ice composition. We intend 
to explore these experimental results in future models. 

For the X-ray desorption rates, we follow the same for- 
mulation as for photodeso rption, cover e d in t he theory 
of .Leger et al.. (,1985l) and iNaiita et al.l (|2001[ ). At this 
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point, it is worth noting that X-ray desorption is the 
least theoretically or experimentally constrained of all 
the non-thermal desorption mechanisms considered here. 
The overall X-ray desorption rate, fc^"", is given by 

f,xr ^ FxRYxBPabsOdXd S~\ (8) 

where, Fxr is the X-ray photon flux in units of pho- 
tons cm~^ s~^, YxR is the desorption yield in units 
of molecules photon"^ and the product, PabsO'd, is the 
effective cross section with Pabs, the probability of X- 
ray absorption by the dust grain. The X-ray desorption 
rate for each individual species, fcf , is calculated accord- 
ing to the fractional abundance of species i on the dust 
grains following Equation ([T]). Her e, we adopt a value 
YxR — 200 from the investigations of lNaiita et al.l ()2001f ) 
and for the e ffective grain cross sect ion we use values from 
the work of iDwek &: SmithI (|1996f ) regarding energy de- 
position into grains by energetic photons i n the energy 
range 10 eV to 1 MeV. iNajita et all (|2001[ ) consider X- 
ray desorption from grains of various compositions and 
morphologies and conclude that both have a significant 
influence on the X-ray desorption yields calculating val- 
ues for YxR ranging between 10 and w 4000 molecules 
photon"^. In this work, we adopt a conservative estimate 
of the yield of 200 molecules photon" -'^ as this is the value 
for the dust morphology which most closely matches our 
simple dust-grain model. Given t he large X-ray lumi- 
nosit ies of T Tauri stars (see e.g. iKastner et al I (|1997l 
12002( 1 )■ we plan a more thorough study on the effects 
of X-ray desorption in protoplanetary disks taking into 
consideration the X-ray energy spectrum as a function of 
disk radius and height and investig ating the full param - 
eter space considered in the work of lNaiita et al.l ()2001[ ). 

2.2.4. Grain-Surface Chemistry 

We use t he grain-surface network fro m lHasegawa et all 
(|1992l l and lHasegawa fc Herbsd (|1993f ) which has 221 re- 
actions involving an additional 9 surface species which 
do not have a gas-phase equivalent (e.g. CH3O). To cal- 
culate the reaction rate coeffici ents, we use the theory 
outlined in detail in ,Hasegawa et al.l (119921 ). The rate 
coefficient for a grain-surface reaction between species i 
and j can be defined as 

kij = Kij {Rdiffii) + RdiffU)) (l/"-d) cm^ (9) 

Here, Kij is the probability that the reaction happens 
upon encounter and is equal to unity for an exothermic 
reaction without an energy barrier. For reactions with an 
activation energy, Ea, and at least one light reactant i.e. 
H or H2, Kij = exp{—2b/ h^/2flEA) where b is the barrier 
thickness and /i = mirrij/ (mi + nij) is the reduced mass 
of the reaction system. This expression is the exponential 
part of the quantum mechanical probability for tunneling 
through a rectangular barrier of thickness, 6. The term, 
Rdi // , is the diffusion rate of an adsorbed species and is 
the inverse of the diffusion time, tdiff, defined as tdiff = 
Ngthop s, where Ns is the total number of surface sites 
per dust grain and t^op is the timescale for an adsorbed 
species to 'hop' from one surface site to another. The 
expression for thop depends on the mass of the species 



and is given by 

_jM^^'ex_p(f^2^^{;E;J^)')s H/H2 

thop — \ \ — 1 f E 

I z^o(0 *^^P ( ~Vr I ^ other species 

(10) 

where Eh{i) « 0.3Ed{i) is the energy barrier between 
surface sites. All other parameters have been defined 
previously. 

3. RESULTS 

We calculate the chemical abundances in the disk as 
a function of disk radius, height and time. The re- 
sults displayed here are extracted at a time of 10^ yr, 
the typical age of visible T Tauri stars with accompa- 
nying protoplanetary disks. Throughout this section, 
fractional abundance refers to the abundance of each 
species with respect to total particle number density. In 
Section 13. 1[ we display and discuss results from model 
PH-I-CRH only, to illustrate the global chemical struc- 
ture and stratification in the disk. Table[2]lists the names 
and ingredients of each model for which we present re- 
sults. Our 'fiducial' model is model PH-f CRH since most 
current chemical models of protoplanetary disks include 
photodcsorption and cosmic-ray induced desorption by 
default. In model CRH we remove photodcsorption to 
investigate the infiuence of cosmic-ray induced desorp- 
tion, in model XD we look at the effects of X-ray des- 
orption only, and in model GR, we investigate the ad- 
dition of grain-surface chemistry to our fiducial model. 
Of course, there are many more permutations of the in- 
gredients which are worthwhile considering in the future 
e.g. X-ray desorption plus grain-surface chemistry. 

3.1. Chemical Structure 

Figure [2] displays the fractional abundances of those 
molecules observed in disks (CO, HCO+, HCN, CN, CS, 
C2H, H2CO and N2H+) as a function of disk radius and 
height, up to maximum radii of 10 AU (left column) and 
305 AU (right column). 

The global abundance distribution of molecules is gov- 
erned by the binding energy of each molecule to dust 
grains and the UV radiation field strength. We see most 
molecules existing predominantly in a molecular layer 
of varying thicknesses at a height, zjr k, 0.3 to 0.5 in 
the outer disk and sa 0.2 to 0.3 in the inner disk with 
freeze out causing depletion in the midplane and pho- 
tolysis causing destruction in the upper layers. CO is 
an exception to this and is abundant (x(CO) « 10"^) 
throughout the majority of the depth of the outer disk 
(> 50 AU) with depletion due to freeze out in the disk 
midplane only occurring beyond a radius of 250 AU. 
In the inner disk (r < 50 AU), gas-phase CO is abun- 
dant in the disk midplane due to its low binding energy 
to the dust grains. In this region, however, we see most 
molecules confined to the 'molecular layer'. In the outer 
disk HCO+ has a peak fractional abundance of ~ lO"^'^ 
to ^ 10"^ throughout most of the disk, mirroring the 
distribution of CO. Within « 50 AU, it is confined to a 
thin layer at a height zjr « 0.3 with x(HCO+) 10"^ 
which coincides with the transition zone where the gas 
composition changes from molecular to atomic hydrogen. 

Gas-phase HCN has a peak fractional abundance of 
x(HCN) ~ 10"'' existing in the molecular layer through- 
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Table 2 

Chemical Models 



CRH PH+CRH XD GR 

Thermal desorption / / / / / 

Cosmic-ray induced desorption / / / 

Photodesorption / / 

X-ray desorption / 

Grain-surface chemistry / 



out the disk. HCN can remain frozen out onto dust grains 
within radii ~ 1 AU of the parent star demonstrating the 
effects of the vastly different desorption energies of CO 
and HCN (960 K and 4170 K, respectively). The distri- 
bution of CN is complementary to that of HCN as it is 
predominantly formed via the photodissociation of the 
latter molecule. Hence, throughout the disk, CN exists 
in a layer above that of HCN with a fractional abundance 
~ 10~^. In the outer disk, CN can survive in the surface 
region, however, the increasing UV field strength in the 
inner disk means that CN is also destroyed by photodis- 
sociation in the disk surface. 

The distribution of the radicals, CS and C2H are sim- 
ilar to that of CN since both are formed predominantly 
via the UV photolysis of larger precursor molecules (e.g. 
H2CS and C2H2). H2CO reaches its maximum frac- 
tional abundance (x(H2C0) ~ 10~®) in the outer disk, 
although, within 10 AU, this value is reduced to ^ 10^^" 
and H2CO is confined to the molecular layer. H2CO is 
returned to the gas phase in the disk midplane within 
« 1 AU. 

The fractional abundance distribution for N2H+ is dif- 
ferent to any of the molecules considered thus far. N2H+ 
reaches its maximum fractional abundance of « 10~^'^ in 
the outer disk only and is present where gas-phase CO 
is depleted e.g. in the disk midplane beyond a radius of 
250 AU. In dense regions, the main destruction mecha- 
nism of N2H+ is via reaction with CO. In the upper lay- 
ers, N2H+ increases in abundance due to the increased 
abundance of both N2 and cations, such as, H3+. Within 
10 AU, x(N2H+) remains less than « lO'^^. 

Figure [3] displays the fractional abundances of those 
additional molecules observed at infrared wavelengths, 
H2O (top), OH (second), CO2 (third) and C2H2 (bot- 
tom), as a function of disk radius and height up to a max- 
imum radius of 10 AU. We display results from within 
10 AU only as infrared emission originates from the inner 
hot, dense disk material. 

Gas-phase H2O is confined to the molecular layer with 
a fractional abundance x(H20) ^ 10^"*. As above, freeze 
out is responsible for depletion in the midplane and pho- 
tolysis for depletion in the upper layers. H2O is returned 
to the gas phase in the midplane at a radius sa 1 to 
2 AU. The distribution of OH is complimentary to that 
of H2O residing throughout the disk in a layer above that 
of gas-phase H2O, reaching a peak fractional abundance 
of ^ lO^"'. The distribution of CO2 is similar to that of 
CO in the inner disk, existing only in the midplane with a 
maximum value of lO"** within a few AU of the star. The 
snow-line for CO2, however, is at the much smaller radius 
of ~ 10 AU (as opposed to w 250 AU). Acetylene, C2H2, 
reaches a peak fractional abundance of ^ 10"*, in the 
molecular layer. C2H2 and similar molecules are formed 
in hotter regions where oxygen is depleted from the gas 



phase via freeze out of oxygen containing molecules onto 
dust grains driving a carbon chemistry and hydrocarbon 
synthesis. 

In Figure O (online only) we display the fractional 
abundances of the gas-phase molecules discussed above, 
along with constituent atoms and grain-surface ana- 
logues where applicable, as a function of disk height at 
radii r = 0.1 AU, 1 AU, 10 AU and 100 AU. 

3.2. Effects of Non-thermal Desorption 

In this section, we discuss the effects of each of our 
non-thermal desorption mechanisms, cosmic-ray induced 
desorption, photodesorption and X-ray desorption, on 
the disk chemical structure (models CRH, PH+CRH and 
XD in Tabled! respectively). We call our control model, 
which includes thermal desorption only, model 0. First, 
we show in Figure |4] those regions of the disk in which 
molecules are depleted if we take into account thermal 
desorption only presenting, as an example, the gas-phase 
CO fractional abundance as a function of disk radius and 
height from model 0. We discuss how efficiently each 
non-thermal desorption mechanism works against deple- 
tion throughout the disk in Sections 13.2.11 to 13.2.31 We 
display the fractional abundances of several gas-phase 
species as a function of disk radius and height comparing 
results from each of our non-thermal desorption models 
in Figure [TU] (online only). 

In Figure |1] there are three notable areas where CO 
is depleted from the gas phase, (i) in the midplane be- 
yond a radius of « 250 AU, (ii) in a layer at a height of 
z/r fa 0.3 and (iii) in the disk surface between a radius of 
a few AU to « 50 AU. In region (i) , the depletion is due 
to freeze out of CO onto dust grains as the dust temper- 
ature here is below the desorption temperature of CO. 
The absence of any non-thermal desorption means that 
CO is completely removed from the gas phase. For region 
(n), the depletion is due to the destruction of CO by UV 
radiation. At this point, the UV field is strong enough to 
dissociate CO into its constituent atoms (atomic carbon 
and oxygen), however, the dust temperature is also low 
enough for the freeze out of molecules which have binding 
energies larger than that of CO e.g. H2O. By a time of 
10^ years, the time at which we extract our abundances, 
atomic carbon and oxygen are trapped in molecules con- 
tained in the icy mantle. Above this height, the dust 
temperature becomes high enough for many molecules to 
desorb thermally thus replenishing the stock of C and O 
to reform CO. In region (iii), a similar effect to that in re- 
gion (ii) occurs due to the decoupling of the dust and gas 
temperatures in the disk surface with the dust temper- 
ature up to two orders of magnitude lower than the gas 
temperature. The thermal desorption rate is dependent 
only on the dust temperature (cx exp{—Ed/Td)) whereas 
the accretion rate depends both on the gas temperature 
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Figure 2. (Continued.) 
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Figure 3. Fractional abundances of H2O (top left), OH (top riglit), 
radius and lieiglit up to a radius of 10 AU. 
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Figure 4. Fractional abundance of gas-phase CO as a function of 
disk radius and heiglit using results from model 0. 



and the number density of dust grains (see Equations [2] 
and [3|) . Since the gas temperature can be much higher 
than the dust temperature in this region, the accretion 
rate can supersede that of thermal desorption so that by 
10^ years, molecules which can survive the intense UV 
radiation field are able to freeze out onto dust grains. 
Specifically, H2O molecules are able to survive as ice on 
dust grains thereby depleting the gas phase of oxygen- 
bearing molecules such as CO. 

3.2.1. Cosmic-ray Induced Desorption 

The left columns of Figures [S] and H] display the frac- 
tional abundances of several molecules and molecular 
ions as a function of disk height at radii r = 10 AU and 
305 AU, respectively, comparing the results from model 




Radius (AU) 

CO2 (bottom left) and C2H2 (bottom right) as a function of disk 

(solid lines) with model CRH (dotted lines). Cosmic- 
ray induced desorption has the smallest effect on the 
gas-phase abundances and only in the outer disk mid- 
plane. The fractional abundances of gas-phase CO and 
H2O are enhanced in the disk midplane at 305 AU in the 
results for model CRH, compared with those for model 
0, although the values reached remain orders of magni- 
tude smaller than those in the upper disk layers. The 
fractional abundance of N2H^ is also enhanced, due in 
part to the release of N2 from dust grains in this region 
by cosmic-rays and, in fact, reaches its maximum frac- 
tional abundance in the disk midplane. As our model 
disk is truncated at 305 AU, we would expect the eff'ects 
of cosmic-ray induced desorption to continue beyond this 
radius. 

3.2.2. Photodesorption 

Photodesorption is the most experimentally con- 
strained non-thermal desorption mechanism we have in- 
cluded in our model. Photodesorption has an effect in 
the upper disk layers where the UV radiation field has 
reached an appreciable strength yet the temperature re- 
mains low enough for freeze out to occur. We compare 
the results from model (solid lines) with those from 
model PH+CRH (dotted lines) in the middle plots of 
Figures [5] and [H Photodesorption enhances the abun- 
dance of molecules in the molecular region of the disk, 
counteracting the depletion effects discussed at the be- 
ginning of Section [XH Clearly seen at a radius of 10 AU 
is the smoothing of molecular abundances throughout 
the middle region of the disk with both model results 
converging higher in the surface. At a radius of 305 AU, 
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there is a smoothing out of abundances counteracting de- 
pletion but also a slight difference in the distribution of 
molecules in the surface regions. The fractional abun- 
dances of CO and CS are enhanced in model PH+CRH 
relative to model 0, whereas those of HCO+, H2O, HCN 
and N2H+ are reduced. This is due to the alteration of 
gas-phase chemistry when photodesorption is included as 
a significant amount of all molecules can remain in the 
gas phase in the mid and upper layers of the disk. Thus, 
those molecules which ordinarily would be frozen out in 
the absence of photodesorption are available to take part 
in gas-phase reactions e.g. N2H+ is destroyed via reac- 
tion with gas-phase CO so that an enhancement in the 
abundance of the latter leads to a corresponding drop in 
that of the former. 

3.2.3. X-ray Desorption 

X-ray desorption is the least theoretically or experi- 
mentally constrained non-thermal desorption mechanism 
we considered, hence, we have used conservative esti- 
mates of molecular yields and thus X-ray desorption 
rates. The right-hand plots in Figures [S] and H] suggest 
that, even using conservative estimates, X-ray desorption 
has the largest effect on gas-phase molecular abundances. 
X-rays, with their higher energy, can penetrate deeper 
into the disk material than UV photons, hence. X-rays 
have an effect in the disk midplane as well as in the molec- 
ular region. At r = 10 AU (top plot), the abundances 
of both H2O and CS are enhanced in model XD rela- 
tive to the results from models 0, CRH and PH+CRH. 
X-ray desorption also acts to smooth out abundances in 
the upper disk, similar to the effects of photodesorption. 
At 305 AU, the effects of the inclusion of X-ray desorp- 
tion are most apparent. The fractional abundances of 
all molecules considered here, with the exception of CS, 
are enhanced in the midplane of the disk, to values com- 
parable with those found in the upper disk regions. In 
fact, it appears that the inclusion of X-ray desorption 
acts to smooth or homogenise the fractional abundances 
of gas-phase CO, H2O and HCO+ throughout the depth 
of the disk. Again, as seen in the results for photodes- 
orption, the gas-phase distributions are altered in model 
XD compared with model due to the alteration of the 
gas-phase chemistry. In the upper disk, the results for 
models PH-t-CRH and XD are similar with the exception 
of the fractional abundance of CS which is enhanced in 
abundance between height of « 70 AU and w 150 AU in 
model XD relative to model PH-I-CRH. 

3.3. Effects of Grain-surface Chemistry 

The addition of grain surface chemistry is expected to 
aid the synthesis of complex organic molecules in regions 
of the disk where significant freeze out has occurred. In 
this discussion, we look at the abundances of small or- 
ganic (saturated) molecules in the outer disk in partic- 
ular. In model GR (see Table [2|), in addition to grain- 
surface chemistry, we also add cosmic-ray induced des- 
orption and photodesorption. Figure [7] shows the frac- 
tional abundances of several small organic molecules as 
a function of disk height at radii, r = 100 AU (left) 
and 305 AU (right), for model PH-hCRH (solid fines) 
compared with model GR (dotted lines). At 100 AU, 
tfie fractional abundance of HCOOH (formic acid) and 



H2CO (formaldefiyde) are enhanced in the disk mid- 
plane in model GR, relative to model PH-I-CRH, with the 
abundance of HCOOH also enhanced in the upper disk 
layer. Most of the organic molecules considered reach 
their peak fractional abundance between 25 and 40 AU. 
Here, the fractional abundances of CH3OH (methanol), 
HCOOCH3 (methyl formate) and CH3OCH3 (dimethyl 
ether) are all enhanced to a value ^ 10"^'^ in model GR, 
orders of magnitude larger than the respective values 
from model PH+CRH. At the very outer edge of our 
disk model, r — 305 AU, the fractional abundances of all 
molecules are enhanced in model GR relative to model 
PH+CRH. Of note is the extreme enhancement seen in 
the abundances of CH3OH, HCOOCH3 and CH3OCH3, 
again by several orders of magnitude, to fractional abun- 
dances which are potentially observable (~ 10"^"'^ to 
~ 10~^^). By fitting observed line intensities of rota- 
tional transi tions in a s e lectio n of molecules with a simple 
disk model, iThi et ahl ()2004[ ) estimate the column den- 
sity of the relatively complex molecule, H2CO, in sev- 
eral protoplanetary disks as lying between ^ 10^^ and 
10^'^ cm^'^ which in the outer disk translates roughly 
to a fractional abundance of ~ 10~^° as the H2 col- 
umn density here is ~ 10^^ cm^^. Complex organic 
molecules in hot cores and dark clouds are routinely 
observed with fractiona l abun dances > 10^^^ (see e.g. 
iHerbst fc van Dishoeckl ()2009t )). We display the frac- 
tional abundances of the molecules discussed in this sec- 
tion as a function of disk radius and height comparing 
results from model PH+CRH and model GR in FigurelTT] 
(online only). 

3.4. Disk lonisation Fraction 

The ionisation fraction in protoplanetary disks is an 
important parameter as it is thought that this drives the 
accretion flow in the disk through the coupling of the 
gas with the strong magnetic fields generated by the sys- 
tem. The required turbulence is generated via magneto- 
rotational instabifities or MRI (Balbus & Hawlcy 1991!; 
iHawlev Sz Balbu^llQQlt ) . For effective accretion, the ion- 
isation fraction is required to exceed a critical level which 
is d ependent on the nature of the star-disk system (see 
e.g. Ilkner fc Nelsm] (l2006h V Re gions in which the ion- 
isation fraction falls below this critical value and where 
effectively magneto-hydrodynamic accretion is switched 
off, are termed 'dead zones'. 

Figure [8] shows the electron fractional abundance as 
a function of disk radius and height within a radius of 
10 AU (left) and 305 AU (right) using the results from 
model PH+CRH. The ionisation fraction varies between 
a minimum value of '-^ 10^ ^■^ in the densest region of the 
disk up to a value of ~ 0.1 in the hottest, most irradi- 
ated surface region closest to the star. We find that we 
attain similar electron abundances throughout the disk 
regardless of chemical model. The ionisation threshold 
required for effecti ve accretion is r elated to the magnetic 
Reynold's number (iGammidliggS *) which must be deter- 
mined in advance of addressing the location of any 'dead 
zones' in our particular star-disk system. The question 
of whether accretion is suppressed in our disk model will 
be considered in detail in a subsequent publication in 
which we also investigate the effects of the recalculation 
of photo-rates and direct X-ray ionisation on the disk 
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Figure 5. Fractional abundances of several gas-phase molecules and molecular ions as a function of disk height at a radius, r = 10 AU 
comparing results from model (solid lines) with each non-thermal desorption model (dotted lines), CRH (left), PH-I-CRH (middle) and 
XD (right). Note that the results for N2H+ from model are too small to appear on our plot. 
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Figure 6. Fractional abundances of several gas-phase molecules and molecular ions as a function of disk height at a radius, r = 305 AU 
comparing results from model (solid lines) with each non-thermal desorption model (dotted lines), CRH (left), PH-I-CRH (middle) and 
XD (right). 



chemical structure and ionisation fraction. 



3.5. Radial Column Densities 

The column density, iV^, at each radius, r, for each 
species, i, is calculated by integrating the number density 
over the depth of the disk i.e. 



ni{r, z) dz 



(11) 



The radial column densities provide an excellent means 
to trace the radial mass distribution in the disk and also 
to compare directly results from each of our chemical 
models to determine the species sensitive to each chem- 
ical process. In Table [3] -we list the column densities of 
various important molecules at radii of 1 AU, 10 AU, 
100 AU and 305 AU for each chemical model and -we dis- 
play the column densities of many of the molecules dis- 
cussed thus far, as a function of radius, up to maximum 
radii of 10 AU (left) and 305 AU (right) in Figure [E] 
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Figure 7. Fractional abundances of several small organic molecules as a function of disk height at radii, r = 100 AU (left) and 305 AU 
(right) for model PH+CRH (solid lines) and model GR (dotted lines). 
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Figure 8. Fractional abundance of electrons as a function of disk radius and height up to maximum radii of 10 AU (left) and 305 AU 
(right). 
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(online only). 

At a radius of 1 AU, the column densities are relatively 
insensitive to the choice of chemical model, unsurprising 
given that most molecules are in the gas phase at this 
radius and the chemistry is dominated in the midplane 
by neutral-neutral reactions and in the surface by photo- 
chemistry. 

At 10 AU, the effects of the inclusion of non-thermal 
desorption become apparent. The column densities for 
those models which include photodesorption are consis- 
tently higher for all molecules with the column densities 
of HCN, CN and CS particularly sensitive. X-ray des- 
orption has a similar effect although we also see a dra- 
matic increase in the column densities of H2O and CO2 
due to the penetrative power of X-rays in this region. 
This is due to the relatively strong X-ray field at this ra- 
dius coupled with the low column density of intervening 
absorbing material (from the disk surface to the mid- 
plane) . Grain-surface chemistry has only a mild effect at 
10 AU on the column densities of the listed molecules. 
For the molecular ions, HCO+ and N2H+, the addition 
of photodesorption causes a rise in the column densities 
of both molecules while the addition of X-ray desorption 
produces a fall in the column density of the former and 
a rise in that of the latter. 

At 100 AU, we see some of the same behaviour as at 
10 AU with photodesorption increasing the column den- 
sities of CO, HCN, CN, CS, C2H, H2CO, H2O, CO2 and 
C2H2 with HCN, CS and CO2 particularly affected. X- 
ray desorption further enhances the column densities of 
these species. Here, we begin to see the effects of grain- 
surface chemistry with HCN, H2O and CO2 casualties of 
the increased synthesis of more complex species. 

Finally, in the outer disk, the effects of cosmic-ray in- 
duced desorption on the column density of N2H+ become 
apparent, so that this molecule is a potential observable 
tracer of this desorption mechanism. We can also see 
how both photodesorption and X-ray desorption act to 
counteract the depletion of CO onto dust grains in the 
outer cold disk midplane thus enhancing its overall col- 
umn density. In the outer disk, observable in the column 
densities is the detrimental effect that the addition of 
grain-surface chemistry has on the column densities of 
HCN and CS in particular, as N and S atoms are in- 
corporated into larger, more complex species, via grain- 
surface reactions. We can also see the dramatic effect on 
methanol due to the inclusion of grain-surface chemistry 
with its column density enhanced by around three orders 
of magnitude. 

3.6. Comparison with Other Models 

A direct comparison with other chemical models of 
protoplanetary disks is difhcult as no two models are 
identical in either their physical basis or their chemical 
networks. Given the plethora of preceding work in this 
field, we limit this short discussion to more recent models 
which are comparable to ours in chemical complexity. 

The work presented here builds upon previous in- 
vestigations into the importance of cosmic-ray in- 
duced desorption and phot odesorpt i on in disks by 
iWillacv fc Lmgol pOOOl ) an dlWillacvl (l200l . the latter 
of which uses the model of iD'Alessio eTaLl (|200l for 
their physical framework. We obtain encouragingly sim- 
ilar results to iWillacv, (,2007. ) in particular, although, her 



Table 3 

Column Densities 



Species CRH PH-I-CRH XD GR 



1 AU 



Ha 


1.9 


25) 


1.9(25) 


1.9(25) 


1.9(25) 


1.9 


25) 


CO 


2.3 


21) 


2.3(21) 


2.3(21) 


2.3(21) 


2.3 


21) 


HCO+ 


1.7 


14) 


1.7(14) 


1.7(14) 


1.7(14) 


1.7 


14) 


HCN 


2.0 


17) 


2.0(17) 


2.0(17) 


2.0(17) 


2.0 


17) 


CN 


2.7 


14) 


2.7(14) 


2.7(14) 


2.7(14) 


2.7 


14) 


CS 


4.5 


12) 


4.5(12) 


4.5(12) 


4.5(12) 


4.5 


12) 


C2H 


4.1 


14) 


4.1(14) 


4.1(14) 


4.1(14) 


4.1 


14) 


H2C0 


1.3 


12) 


1.3(12) 


1.3(12) 


1.3(12) 


1.3 


12) 


N2H+ 


5.5 


10) 


5.5(10) 


5.3(10) 


5.2(10) 


5.3 


10) 


OH 


1.7 


16) 


1.7(16) 


1.7(16) 


1.7(16) 


1.7 


16) 


H2O 


1.7 


21) 


1.7(21) 


1.7(21) 


1.6(21) 


1.7 


21) 


CO2 


4.6 


20) 


4.6(20) 


4.6(20) 


4.6(20) 


4.6 


20) 


C2H2 


1.8 


15) 


1.8(15) 


1.8(15) 


1.8(15) 
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primary objective was the investigation into deuterated 
species in disks. As such, the reaction network of the non- 
deuterated species was truncated to accommodate the 
additional reactions involving deuterium and deuterium- 
containing molecules. Our work also differs in that our 
physical model includes X-ray heating and also explic- 
itly determines the gas temperature which can decouple 
from the dust-grain temperature in re gions where cooling 
via ga s-grain collisions is inefficient. iWillacv fc Woodsl 
(pOOa ) also use their deuterated reaction network and 
adapt their model to investigate the chemical structure 
of a disk within 30 AU of the central star. Again, we 
achieve similar results although differences in the set of 
molecular desorption energies used manifests as differ- 
ences in the positions of 'snow lines' for molecules such 
as HCN. Also differing prescriptions for the UV radiation 
field in the disk leads to different distributions of radicals 
such as C2H and CN. 

Photodesorption and cosmic-ray induced desorption 
are now routinely included in m odern chemic a l mod - 
els of protopl a netary disks (e.g. iWoitke et al.l ()2009D : 
IHenning et all (pOm) ). X -ray desorption and grain- 
surface chemistry have both been i ncluded in work by by 
other groups (e.g. iSemenov et al.l (|2008[ ) : iHenning et al.l 
(|2010D ') although not explicitly investigated given the 
theoretical uncertainty behind the exact mechanism of 
X-ray desorption for the former process and the usual 
truncation of chemical networks for the latter. Our work 
presented here and subsequent follow-up publications on 
X-ray desorption and grain-surface chemistry, we hope, 
will go some way to addressing this. 

4. SUMMARY 

In this work, we have presented a selection of results 
from our high-resolution combined chemical and physi- 
cal model of a protoplanetary disk surrounding a typical 
T Tauri star, constructed in order to trace the physical 
and thus chemical structure on small scales. We use a 
protoplanetary disk model in which the gas density and 
temperature distributions are obtained self-consistently 
and UV and X-ray irradiation by the central star is calcu- 
lated by solving the radiative transfer equation. To this 
we applied a large comprehensive chemical network in- 
cluding gas-phase chemistry, gas-grain interactions and 
grain-surface chemistry. We investigated the effects of 
each non-thermal desorption mechanism thought to be 
important in disks: cosmic-ray induced desorption, pho- 
todesorption and X-ray desorption. We also added a 
large grain-surface network to investigate the effective- 
ness of grain-surface reactions on the synthesis of rela- 
tively complex organic molecules. 

Using the results from model PH-I-CRH (see Table H]) 
extracted at a time of 10® years, we find that the disk 
chemical structure closely mirrors the disk physical struc- 
ture with the freeze out of molecules onto dust grains 
creating an icy mantle in the cold, dense midplane and 
an abundance of molecules in a layer above the mid- 
plane created through sublimation and the resulting rich 
gas-phase chemistry. In the disk surface, the molecular 
abundances drop as the UV and X-ray radiation fields 
peak in strength dissociating molecules and ionising both 
molecules and atoms. The resulting disk ionisation frac- 
tion increases with increasing disk height. There is simi- 
lar stratification in the radial direction as the increasing 



temperature drives the evaporation of molecules. The 
temperature dependence of the desorption energies re- 
sults in a unique 'snow-line' for each molecules with more 
volatile molecules returned to the gas in the midplane at 
larger radii than less volatile ones. In particular, both 
HCN and H2O remain frozen onto dust grains to within 
« 1 to 2 AU of the central star. 

The addition of cosmic-ray induced desorption has only 
a small effect on the gas-phase abundances in the outer 
disk midplane although since our model is truncated at 
305 AU we would expect that this effect will continue at 
radii beyond this value. Photodesorption, the most ex- 
perimentally constrained of the mechanisms considered 
here, has a larger effect although only in the molecular 
and surface regions of the disk where there is an appre- 
ciable UV fiux. It is especially effective at enhancing the 
gas-phase abundances of non-volatile molecules such as 
II2O. X-ray desorption, has the largest effect, smooth- 
ing the abundances of gas-phase species throughout the 
height of the disk and acting to homogenise the fractional 
abundances. However, X-ray desorption is the least the- 
oretically or experimentally constrained and thus our re- 
sults must be treated with caution, pending further in- 
vestigation. 

The addition of grain-surface chemistry also yields 
some encouraging results worthy of revisiting and fur- 
ther study. In the outer disk, where the freeze out of 
molecules is most prevalent, the abundances of relatively 
complex organic molecules e.g. CH3OH, HCOOCH3 and 
CH3OCH3 are enhanced to potentially observable values 
when grain-surface chemistry and photodesorption are 
included in our model. Thus, the observation of rota- 
tional transitions of these, and related, species in pro- 
toplanetary disks should provide an excellent means of 
testing grain-surface chemistry theory. Due to limita- 
tions in existing facilities, the most complex molecule 
observed, as yet, in disks is H2CO. ALMA, however, will 
have the sensitivity and spectral resolution necessary to 
observe rotational transitions in these minor species. In- 
deed, preliminary synthetic spectra we have calculated 
suggest that the rotational transition lines of methanol 
are enhanced above the detection threshold of ALMA 
when grain-surface chemistry is included in our model 
(Walsh et al., in preparation). 

We have shown that running models of this nature, 
in which we test experimental data and theory, varying 
the different chemical ingredients, are necessary for dis- 
entangling the different physical infiuences on the chem- 
ical content of protoplanetary disks. The influence of 
the physical conditions and processes on the molecular 
content is also a strong function of radius and so high- 
resolution models, which trace the chemical structure on 
small scales, are also preferred. In this brief overview of 
our model, we have shown that X-ray desorption and 
grain-surface chemistry can have a powerful effect on 
the molecular content of disks and we intend to expand 
upon the work presented here with follow-up papers on 
both chemical processes. Although in this paper we have 
shown that the distribution and radial column densities 
of particular molecules are sensitive to the inclusion or 
omission of certain chemical processes, in order to test 
the viability of using these molecules as tracers we must 
compute the radiative transfer in the disk and directly 
compare our results with observations. This work will be 
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reported in a subsequent paper (Walsh et al., in prepa- 
ration). 
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Figure 9. Fractional abundances as a function of disk height at radii, r = 0.1 AU (top), 1 AU (second), 10 AU (third) and 100 AU 
(bottom). Grain-surface (ice) abundances are represented by dotted lines. Figure available online only. 
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Figure 10. Fractional abundances as a function of dislc radius and height for several gas-phase molecules for model CRH (left), PH+CRH 
(middle) and XD (right). Figure available online only. 
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Figure 11. Fractional abundances as a function of disk radius and height for several organic molecules for model PH+CRH (left) 
GR, (right). We display results from model GR only for HCOOCHa and CH3OCH3 (bottom row) as model PH+CRH produces neglij 
fractional abundances (< 10~^^) for both of these molecules throughout the disk. Figure available online only. 
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Figure 12. Radial column densities of various gas-pliase molecules up to maximum radii of 10 AU (left) and 305 AU (right) for each 
chemical model. Grain-surface (ice) column densities are represented by dotted lines. Note that results from models and CRH (red 
and green lines, respectively) are very similar hence for many molecules these lines are virtually indistiguishable. Note also for CO2, in 
particular, the results from models 0, CRH and PH+CRH within 10 AU arc almost identical. Figure available online only. 
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APPENDIX 
DISK PHYSICAL STRUCTURE 
Density 

The number density (cm~'^), as a function of disk radius and height, is displayed in the top left and right panels 
of Figure [1] (within a radius of 10 AU and 305 AU, respectively) and, as a function of z/r (i.e. disk height scaled by 
radius) at a range of radii, in the top left panel of Figure [T31 The density decreases as a function of increasing disk 
radius and height with the densest region of the disk found in the disk midplane close to the star 10^^ cm"'^) and 
the most diffuse, in the disk surface at large radii (~ 10^ cm~'^) so that the density range in our model covers almost 
10 orders of magnitude. The density gradient in the vertical direction becomes more extreme with decreasing radius. 
At a radius of 0.1 AU, the density ranges from ^ 10^ cm""^ in the disk surface to ^ 10^^ cm~^ in the midplane, a 
change of seven orders of magnitude over a distance of w 0.03 AU. In contrast, at a radius of 305 AU, the density 
varies by around a factor of 100 between the midplane and the surface, over a distance of « 250 AU. 

Gas and Dust Temperature 

Corresponding graphs for the gas temperature are shown in the bottom left and right panels of Figure [T] and the 
bottom left panel of Figure [T51 The gas temperature increases as a function of increasing height and decreasing radius 
with the hottest region found in the disk surface close to the star (~ lO"' K), and the coldest found in the disk midplane 
in the outer disk (~ 10 K). Within several AU of the star, the temperature increases near the disk midplane due to the 
influence of viscous heating. In the disk surface, the gas and dust temperatures decouple with the dust temperature 
the lower of the two values (represented by dotted lines in Figure fT3|) . At low densities, gas-grain collisions become 
ineffective so the gas cools via radiative line transitions. At a radius of 10 AU, the difference in the gas and dust 
temperatures in the disk surface is more than an order of magnitude at « 2000 K and « 100 K, respectively. The 
vertical temperature gradients in the disk, although large, are not as extreme as those of the density. 

UV and X-ray Radiation Fields 

The middle panels of Figure [T51 show the UV and X-ray fluxes (top and bottom, respectively) as functions of z/r at a 
range of radii. The UV flux is given in units of the average interstellar radiation field (ISRF) (1.6 x 10~^ erg cm~^ s~^). 
The extinction coefhcient in the disk is very large for UV radiation thus the UV flux has an incredibly steep gradient 
in the vertical direction with the flux in midplane at all radii « 0. In the disk surface, the UV radiation field decreases 
in strength with distance from the star. At a radius of 0.1 AU, the UV flux has a value ~ 10^ times that of the ISRF 
whereas at 305 AU, this has decreased to ^ 10 times the interstellar value. 

The X-ray flux shows a similar behaviour, however. X-ray photons, with energies ~ keV, are more energetic than 
UV photons and so have greater penetrative power. Hence, the X-ray flux has a less steep gradient in the vertical 
direction than the UV flux. In the outermost regions of the disk, r > 10 AU, the X-ray flux is attenuated from a value 
~ 0.1 to ~ 10 erg cm~^ s~^ in the disk surface to ~ 10~^ to ~ 10~^ erg cm~^ s~^ in the disk midplane so that there 
remains a small, yet significant. X-ray flux where the UV flux is essentially zero. 

Cosmic-ray and X-ray lonisation Rates 

Finally, the cosmic-ray and X-ray ionisation rates are displayed in Figure [13] (top and bottom right panels, respec- 
tively) as functions of z/r at a range of radii. Note the linear scale used for the cosmic-ray ionisation rate. Beyond 
r 10 AU, the cosmic-ray ionisation rate is constant (at the interstellar value of ~ 10~^^ s~^) throughout the height 
of the disk as cosmic-ray particles, with energies ~ MeV, are highly penetrative. In the inner disk, e.g. r — 1 AU, the 
column density of material through to the disk midplane is large enough for some attenuation of cosmic-rays, however, 
the midplane ionisation rate at this radius in our model is still relatively high at « 60 % of the interstellar value. At 
even smaller radii, r ~ 0.1 AU, the midplane cosmic-ray ionisation rate is further reduced to « 20 % the rate in the 
disk surface. 

The X-ray ionisation rate, unsurprisingly, mirrors the X-ray flux. In the outer disk, r > 10 AU, the X-ray flux 
is significant enough to result in an X-ray ionisation rate ~ lO^^*' to ~ 10~^^ in the midplane, however, here, 
cosmic-ray ionisation (~ lO"^'' s^^) dominates. The X-ray ionisation rate is significantly higher than the cosmic-ray 
ionisation rate in the middle to surface layers of the disk, with values ranging from ^ 10^^'^ s^^ at the maximum 
radius of our disk model to ~ 10~^ s~^ at the minimum. 
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Figure 13. Number density (top left), temperature (bottom left), UV flux (top middle). X-ray flux (bottom middle), cosmic-ray ionisation 
rate (top right) and total X-ray ionisation rate (bottom right) as a function of disk height at radii, r = 0.1 AU, 1 AU, 10 AU, 100 AU and 
305 AU. The gas and dust temperatures decouple in the disk surface with the dust temperature represented by the dashed lines. The UV 
flux is given in units of the flux of the interstellar radiation field or ISRF (1.6 X 10~^ erg cm~^ s~^). Note that the scale used for the 
cosmic-ray ionisation rate is linear. 



